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Abstract 

Relativistic particles in the Kepler and Coulomb potentials may have trajectories that are qual- 
itatively different from the trajectories found in nonrelativistic mechanics. Spiral scattering tra- 
jectories were pointed out by C. G. Darwin in 1913 in connection with the relativistic Rutherford 
scattering of classical charged particles. Relativistic trajectories are of current interest in connec- 
tion with Cole and Zou's computer simulation of the hydrogen ground state in classical physics. 
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I. INTRODUCTION 



The mechanics of a relativistic particle in a potential is usually treated as an afterthought 
in classical mechanics textbooks. These textbooks usually begin with an extensive treat- 
ment of nonrelativistic mechanics, then the energy and momentum of relativistic particles are 
noted, and the corrections to the orbits of nonrelativistic particles are sometimes mentioned. 1 
However, relativistic particles in the familiar Kepler or Coulomb potential V(r) = —a/r can 
have trajectories that are qualitatively different from the trajectories found from nonrela- 
tivistic mechanics, and these unfamiliar trajectories are not mentioned in the textbooks.— 
For example, a relativistic particle in a 1/r potential can spiral into the force center (while 
conserving mechanical energy and angular momentum). This behavior occurs because a 
small increase in the velocity near the speed of light c can lead to a large increase in the 
mass mj (1 — v 2 /c 2 ) 1 / 2 so that an increase in the kinetic energy will compensate the decrease 
in the potential energy as the radius decreases, thus keeping the total energy constant. In 
addition, the increase in the linear momentum will compensate the decrease in the radius to 
keep the angular momentum constant. The existence of such trajectories has played a role 
in recent research.^ 

The potential V(r) = —a/r (a > 0) appears in the Kepler problem of gravitational 
physics and in the attraction of point charges in electrostatics. The theory of gravitation 
finds its relativistic form in general relativity and electrostatics has a natural extension into 
electrodynamics. However, we will not discuss these extended theories. The problem we 
consider has no curved spacetime and no radiation; it is the relativistic mechanics of a 
particle in a 1/r potential. 

As a first example of the unfamiliar nature of some relativistic particle trajectories, we 
show that if the angular momentum is too small, there are no circular orbits for a relativis- 
tic particle in a 1/r potential. In contrast, in nonrelativistic mechanics, such a situation 
never occurs; there are always nonrelativistic circular orbits unless the angular momentum 
vanishes. 

We follow the traditional procedure for the classification of trajectories for a particle, 
noting the orbits that undergo a qualitative transformation in the nonrelativistic limit and 
those that do not. The equations of the trajectories are then obtained. We find the familiar 
scattering trajectories and the familiar rosette shapes that appear in texts on old quantum 
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theory^ and also C. G. Darwin's spiraling trajectories,— which do not appear in textbooks. 
Finally, we mention the research context in which these spiral trajectories have recently 
become of interest. 



II. THE MECHANICAL PROBLEM 

We compare the nonrelativistic and relativistic behavior of a particle in the potential 
V(r) = —a/r. A particle in such a potential experiences a force F = -VV(r) = —ra/r 2 . 
If we use the nonrelativistic momentum mv for a nonrelativistic particle, Newton's second 
law becomes 



|(mv) = -^f. (1) 



For a relativistic particle, the momentum is mv/ y 1 — v 2 /c 2 , which gives 

d , mv x a A 

If we take the dot product of Eqs. (Q) and (0) with the velocity v and integrate with respect 
to time, we obtain conservation of energy. The nonrelativistic energy £ nr is expressed as 

£nr = \mV 2 - -, (3) 

2 r 

while the relativistic energy E is 

^ 9 m c 2 a 

E = £ + mc = - . 4 

y/1 -v 2 /c 2 r 

Here we have written E = £ + mc 2 so that £ is the energy difference from the particle rest 
energy. Also, because the 1/r potential gives a central force, the angular momentum L is 
conserved, giving in the nonrelativistic case 

L nr = r x (mv), (5) 

and in the relativistic case 

We will use the conservation laws to determine the trajectories. As an introduction to the 
nature of the differences that appear in relativistic mechanics, we first consider the case of 
circular orbits. 
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III. LIMITING ANGULAR MOMENTUM FOR RELATIVISTIC CIRCULAR OR- 
BITS 



For circular orbits, the particle displacement r from the center of the potential is perpen- 
dicular to the velocity v, and the angular momentum L is perpendicular to the plane of the 
orbit. The magnitude of the angular momentum in the nonrelativistic case is 

L nr = mrv, (7) 

and in the relativistic case 

L = (8) 

y/l ~ V 2 /C 2 

For circular orbits in the nonrelativistic case, we have 

m- = -, (9) 



and in the relativistic case 



m v 2 a 



, ,. (10) 

y/l - v 2 /c 2 r r 2 

If we now use the angular momentum to remove either r or d, we find for the nonrelativistic 

case 

a L 2 
v = — or r=-*£. (11) 
L nr ma 

Because both v and r can take all values between zero and infinity in nonrelativistic me- 
chanics, it follows that any value of the angular momentum L nr can lead to a circular orbit. 
The same procedure in the relativistic case leads to 

v = H (12a) 

or 

'= — [l"(f)T- < 12b > 
ma Lc 

In the relativistic case, we have an upper limit for the speed v — c, and accordingly a lower 
limit on the magnitude of the angular momentum for a circular orbit. As seen in Eq. (|12|). 
we must have 6 

L > - c (13) 

Thus there is a qualitative distinction for circular orbits between nonrelativistic and rela- 
tivistic mechanics for the potential V(r) = —a/r, and only for this potential. 7 The limiting 



value of the angular momentum in the relativistic case is L a = a/c. In contrast, in non- 
relativistic mechanics there is no limiting speed and thus no lower limit on the angular 
momentum for circular orbits. 

The actual value of the limiting angular momentum L a for the 1/r potential depends on 
the magnitude of a. For an electron of charge e in the field of a nucleus of atomic number 
Z, a = Ze 2 , and the limiting angular momentum can be written as L a = Zm(e 2 /mc 2 )c = 
Z(e 2 /hc)h. Thus for an electron in the Coulomb potential of hydrogen , the limiting angular 
momentum L a = (e 2 /hc)h ~ (1/137)^,. However the reader should not be distracted by an 
angular momentum that is close to that given by Planck's constant. Our discussion involves 
only relativistic classical mechanics, not quantum mechanics. 



IV. CLASSIFICATION OF TRAJECTORIES 



The classification of orbits for the potential V(r) = —a/r reflects the fact that some rela- 
tivistic trajectories are qualitatively different from those found in nonrelativistic mechanics. 
The classification can be made following the usual procedures of classical mechanics.— The 
motion is confined to a plane so that the velocity can be written in terms of polar coordinates 
as r = fr + 9r9. Then we use the energy expression (4) and remove the 9 dependence in the 
velocity in favor of the angular momentum L to obtain a first-order differential equation in 
the radial variable r as a function of time. We will carry out the procedure for the relativistic 
case and then take the nonrelativistic limit to connect with the familiar Kepler orbits. 

In terms of polar coordinates, the relativistic angular momentum is 

mr 2 9 (M) 



yjl - (r 2 + r 2 9 2 )/c 2 



where we choose the orientation so that L and 9 are positive. If we solve Eq. (fPfj) for 9 and 
substitute the result into Eq. (JU), we find 



mc 2 a 



E = S + mc 2 = ; --. (15) 

y/1 - (f 2 /c 2 ) - L 2 (l - r 2 /c 2 )/(m 2 r 2 c 2 + L 2 ) r 



The solution of Eq. (|15|) for r 2 is 



•2 2 

r = c 



, L 2 mc 2 
' + m 2 r 2 c 2 ^E + (a/ry 



(16) 
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In the following, it will turn out that the limiting angular momentum, L a = a/c, that 
appears in the analysis of circular orbits, represents a crucial limiting value. It is useful 
to have certain constraints on the energy which are associated with this limiting angular 
momentum. Because the particle velocity is less than c, we see from Eqs. (II 4|) and (|15|) that 



mr 2 9 rmc r , a-, , s 

■ < ; = -(E+-). 17 

y/1 ~ V 2 /c 2 y/1 - V 2 /c 2 C K V > 



Equation (|T7|) requires that 

r, a N a Er . . 

L <-(E+-) or L < — . 18 

c r c c 

Because r is positive, it follows that if L > a/c, then E = £ + mc 2 > 0. Thus orbits 
with angular momentum larger than L a must have positive total (relativistic) energy. In 
particular, bound orbits with L > L a (those that do not plunge into the potential center) 
cannot have arbitrarily small values of kinetic plus potential energy; rather the kinetic plus 
potential energy £ must be larger than —mc 2 . 

The classification of the relativistic orbits can be made by noting that the function on the 
right-hand side of Eq. (fT7)|) must satisfy the requirement < r 2 < c 2 . It is straightforward 
to transform this condition to 

-L 2 c 2 < < (E 2 — m 2 c 4 )r 2 + 2Ear + {a 2 - L 2 c 2 ). (19) 

Because all quantities are real, the first inequality in Eq. (|T9*j) holds automatically and is 
not important. The second inequality can be treated by plotting the parabolic function 
Y(r) = (E 2 — m 2 c 4 )r 2 + laEr + (a 2 — L 2 c 2 ) versus r. The allowed orbits correspond to the 
regions of positive r where Y(r) > and the turning points occur when Y(r) = 0: 

Ea ± ^E 2 a 2 + (m 2 c 4 - E 2 )(a 2 - L 2 c 2 ) 

'"turning— point n ~ A " \ ) 

m A & — E z 

In the nonrelativistic limit c —>■ oo with E — mc 2 = £ — > £ nr , the inequality ()19j) becomes 

-L 2 nr /m < < £ nr r 2 + ar - L 2 nr /2m, (21) 

where the first inequality holds for any real L nr . The parabolic function y(r) = £ nr r 2 + ar — 
L 2 ir /2m can be plotted versus r with the allowed orbits corresponding to regions of r > 
and y(r) > 0; the turning points occur when y(r) = 0, 



A relativistic circular orbit corresponds to the square root in Eq. (|2"Uj) equal to zero so 
that both the inner and outer turning points are at the same radius r. The vanishing square 
root gives a connection between E and L for a circular orbit, 



2 / -, / « \ 2 



£ = mc ^l-(_j, (23) 

If we substitute Eq. into Eq. (J2UJ), we obtain a radius r, which is the same result obtained 
in Eq. (jl2b|) . Also, if we take the nonrelativistic limit of Eq. (|23)L we obtain 



£ = E-mc 2 = mc 2 Jl-^—J - mc 2 91 --— = £ nr , (24) 



nr 



which is the nonrelativistic expression for a circular orbit. We note that Eq. ()23|) gives a 
limit E > for a circular orbit, which agrees with the result given below Eq. (JTHJ). However, 
there is no lower bound for the energy in the nonrelativistic approximation appearing in 
Eq. ()24j) when we take the nonrelativistic limit L nr —>■ (and find E nr — > — oo). 

The general character of the turning points in the nonrelativistic equation (}2"2"|) is con- 
trolled by the sign of the single parameter £ nr , with L nr = a unique special case. In the full 
relativistic treatment, the character of the turning points in Eq. (120}) depends on the signs 
of both i?a/(m 2 c 4 — E 2 ) and (m 2 c 4 — E 2 )(a 2 — L 2 c 2 ), where there is a sign change in the 
second term depending on the magnitude of L. When L > L a so that (a 2 — L 2 c 2 ) < 0, the 
turning point analysis is qualitatively the same for both the relativistic and nonrelativistic 
cases. However, if the angular momentum is small but non-zero, < L < a/c, so that 
(a 2 — L 2 c 2 ) is positive, then the turning point analysis allows new possibilities which do not 
appear in the nonrelativistic case. We will see these new trajectories when we obtain the 
orbit equations. 

V. ORBIT EQUATIONS 

The orbit equations r(6) can be found using the traditional analysis. If we write the 
relativistic particle momentum in the form 

- m(rr + Orb) mrr ~L 

p = rp ' + ep ° = 7T^W = 7T^ + 9 7' (25) 

then the relativistic particle energy E can be written as 



I Ui 

E = ^/p 2 c 2 + m 2 c 4 - -, (26) 



or 

T 2 2 

a \2 2 2, 24 22, -^C 



/ tti i \ 22, 24 22,-^^ i 24 / 7\ 

(-EH — J =pc +mc = p r c H — h m c . (27) 

(28) 



r 

Now from Eq. (J2SJ) 

p r r ldr 
pe r6 rd6' 

and 



r 2 d0' 

Thus Eq. ([2*7]) becomes 



Pr = -,-Tn- ( 2 9) 



a, 2 / L dr N 2 9 £ 2 c 2 



r r 2 dd r 2 

If we introduce the usual r-inverse variable, s — 1/r, so that 

(is 1 c?r 

Eq. (J3TIJ) becomes 

-(iSx2 2r2 2 , _2„4 



(31) 



(£ + as) 2 = L 2 c 2 Gr + ^'c 2 + ™ 2 c*. ( 32 ) 



Next we differentiate Eq. ()32j) with respect to and divide through by ds/dO to obtain 

^2 + 



<i 2 s r / a x 2l -Eq; 



dB l 



1 -(^)^-tB = - (33) 



Equation ()33|) is a second-order linear differential equation and can be easily solved. Then 
we substitute the solution into the first-order differential equation (}3*2*|) to determine some 
of the integration constants. 

The connection between position and time can be obtained from Eq. (jlfij) for f = dr/dt, 

r 

f^f-l h -, r mC 2 ]2\-l/2 , 

t = J ^ + riW^E + (a/ry > ' [ ) 

Alternatively, the angular momentum equation can be combined with the orbital equation. 
Thus Eqs. (fTlj) and (fT5j) can be combined to give 

^ L (35) 



c 2 E + a/r' 
so that from 6 = dd/dt, we have 



t= q/o (36) 



where r(#) is to be evaluated from the orbital equations. 



VI. SOLUTION OF THE ORBITAL EQUATIONS 



We obtain different solutions of Eqs. (J3*2|l and (JHHj) depending on the value of the angular 
momentum L. If L > a/c, we have 



1 E 2 L 2 c 2 -m 2 c\L 2 c 2 -a 2 ) r T Ta7% fn „ x1 Ea 

s =- r = v — — - cos [y 1 " fe } ( " " " o)] + w^?- (37) 

If L = a/c, we have 

8 = I = i(£)(,-«.). + =^-^. ,38, 
r 2 v a 2_&a 

For L < a/c, we have 



1 /m 2 c 4 (a 2 - L 2 c 2 ) + £ 2 L 2 c 2 r /To". ~ „ . £a 



• = 7 = V c ° 8h [ V ( Z?) ^ 1( " - - CT^' (39) 

Finally, if L = 0, the orbit is a straight line 6 = 9q. 

In the nonrelativistic limit 1/c — > 0, only the first of these solutions (L > a/c) does 
not degenerate into a straight-line orbit. The relativistic orbits in Eq. (J37j) are the bound 
rosettes and unbounded scattering orbits which usually appear in textbooks.— Thus if we set 
E — £ + mc 2 and take the limit 1/c — > in Eq. ([37)1. we find 



1 ma / 2L 2 .£ nr ma , ln . 

a = r = iiX + ^ ms{e - " ) + ii; (40) 

which is the standard nonrelativistic result for the hyperbolas, parabolas, and ellipses of 
Kepler orbits. 

The various types of orbits are sketched in Figs. 1-5, showing bound and unbound rela- 
tivistic trajectories. The unbound orbits, shown in Figs. 1-3, have a total relativistic energy 
at least as large as the particle rest energy, E = £ + mc 2 > mc 2 . These orbits extend to 
spatial infinity in at least one of their time limits. For L > a/c, these orbits are familiar 
scattering orbits as shown in Fig. 1, where the orbit looks like a parabolic or hyperbolic 
nonrelativistic orbit. However, for L > a/c but close to a/c, we are reminded that relativity 
changes the familiar nonrelativistic scattering orbits. In Fig. 2 we see that the scattering 
trajectory loops around the center of the potential before receding to infinity. These loops do 
not occur for the parabolic or hyperbolic orbits of nonrelativistic scattering. As L decreases 
and becomes closer to a/c, the number of times the orbit loops around the potential center 
increases. When L < a/c, the character of the trajectory changes. The orbital looping now 



continues all the way to the center of the potential; the unbound orbit extends to infinity 
in only one time direction and spirals into the potential center in the other (see Fig. 3). In 
contrast, the parabolic and hyperbolic orbits of nonrelativistic mechanics never reach the 
center of the potential, unless L = 0. 

The bound orbits have total energy E smaller than the particle rest energy, E = S+mc 2 < 
mc 2 , corresponding to negative values of £ = E — mc 2 . In nonrelativistic mechanics, all the 
bound orbits are ellipses unless L nr = 0. For L > a/c, the relativistic orbits take the rosette 
shape of precessing ellipses shown in Fig. 4. These orbits have two turning points. The 
orbits are familiar as the Sommerfeld relativistic orbits of old quantum theory— and reduce 
to the familiar ellipses in the nonrelativistic limit. For L < a/c, the relativistic orbits spiral 
out from the center of the 1/r potential and back into the center at early and late times as 
shown in Fig. 5. In the nonrelativistic limit, these orbits become straight-line trajectories 
for L nr = 0. 

From Eqs. (}34j) and (}36|) . we find that all the trajectories involve finite time changes for 
finite changes in r or 9. For the trajectories with L < a/c (which spiral into the potential 
center), the particles arrive at the potential center in a finite time starting from a finite 
radius. This result can be seen from Eq. (pMj) which connects the change in the radius r 
with the elapsed time. For L = a/c, an expansion for small r gives an integral of the form 



which is well behaved at the lower limit. If L < a/c, an expansion for small r gives 



The square root is real because of the condition L < a/c, and again the integral is finite. 
Thus all the trajectories involve only finite times in the vicinity of the 1/r potential center. 

VII. RELEVANCE OF THE TRAJECTORIES TO CURRENT RESEARCH 

The existence of classical relativistic scattering trajectories that spiral into the center 
of the 1/r potential while preserving energy and angular momentum was pointed out by 
Darwin 2 in 1913. D arwin's analysis was made in connection with electron scattering based 




(41) 



o 




(42) 
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on Rutherford's model of the nuclear atom. In particular, Darwin was concerned with the 
spiral trajectories and suggested that they should not occur in nature because they would 
cause chemical transmutation. He speculated that "There must therefore be some way by 
which the electron can escape from the extreme neighborhood of the nucleus."— He suggested 
that the presence of these unobserved spiral trajectories indicated a further failure of classical 
electromagnetic theory to describe experimental observation. Darwin's work was published 
just before Bohr proposed his stationary-state model of the atom. 

Although Darwin was concerned primarily with the Rutherford scattering problem of 
relativistic mechanics, he estimated the radiative corrections for his scattering trajectories. 
It is interesting that a relativistic treatment of the old classical problem of radiative atomic 
collapse gives a different result from the familiar nonrelativistic treatment. In both the 
nonrelativistic and relativistic analyses, the radiative loss of mechanical energy as a charged 
classical particle moves in a Coulomb potential is such as to make the orbit more circular. In 
the nonrelativistic treatment, the charged particle radiates an infinite amount of energy as 
it spirals into the center at ever higher speed and ever lower energy S nr = —Ze 2 /2r — ► — oo. 
However, in the relativistic case, the radiative loss of energy ends when the particle speed 
in the circular orbit has reached the limit v — > c, the angular momentum in Eq. (fTHjl has 
reached L — > Ze 2 /c, the radius in Eq. (JT2j) has shrunk to r — > 0, and the total energy in 
Eq. has reached E — > 0. Thus a charged particle of rest energy mc 2 that starts at a 
large radius and small velocity with total energy E = mc 2 will radiate away the energy mc 2 , 
and not the divergent radiation energy found in the nonrelativistic treatment. 

Today, physicists are not concerned about Darwin's spiral trajectories nor the finite radi- 
ation energy loss given by a relativistic treatment of atomic collapse. Of course, most physi- 
cists are not aware of their existence, because quantum mechanics has changed our views of 
the atom and has made such awareness unnecessary. However, occasionally some researchers 
are interested to see just how far one can push a classical or semi-classical interpretation 
of atomic physics. It has been suggested^ that the inclusion of classical electromagnetic 
zero-point radiation might allow an understanding of electron motion as a sort of Brownian 
motion which avoids the problem of atomic collapse. Thus an electron might indeed lose 
energy through radiation while accelerating in the electric field of the nucleus, but the elec- 
tron might also pick up energy from the zero-point radiation. The balance between these 
two effects might account for atomic structure. 
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Recently Cole and Zou2> have carried out simulations of the classical hydrogen ground state 
assuming classical zero-point radiation. They conclude that the electron neither falls into the 
nucleus nor is ionized, but rather assumes a stationary probability distribution. This ground 
state probability distribution approximates that given by the Schrodinger equation. This 
agreement is remarkable and requires confirmation and further understanding. However, 
in contrast to these computer simulation results are various attempts to calculate the elec- 
tron Brownian motion within the basic zero-point model using nonrelativistic mechanics.— 
These analytic calculations suggest that the electron would not fall into the nucleus due 
to radiation emission, but rather would be ejected from the atom (self-ionization) due to 
the absorption of excessive zero-point energy in the plunging orbits of very small angular 
momentum. However, as we have emphasized here, the nonrelativistic orbits of very small 
angular momenta represent the region of failure of the nonrelativistic approximation. For 
small angular momentum L < a/c, the correct relativistic orbits spiral in toward the nucleus 
while conserving the total energy and angular momentum. This behavior will lead to a to- 
tally different interaction with zero-point radiation from that assumed in the nonrelativistic 
calculations.— Hence, the validity of the self-ionization conclusions in the analytic work are 
questionable. 4 

We see that nearly a century after Darwin's work, a knowledge of the orbits of a relativistic 
particle in a 1/r potential again seems significant. An understanding of the possibilities and 
limitations of a classical or semi-classical treatment of hydrogen depends upon a correct 
treatment of the mechanical trajectories of a relativistic particle in the Coulomb potential. 
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Figure Captions 

Fig. 1. Sketch of an unbound orbit of total energy E = £ + mc 2 > mc 2 and angular 
momentum L » a/c in a 1/r potential centered on the origin. The orbit is similar to a 
parabolic or hyperbolic nonrelativistic scattering orbit. 

Fig. 2. Sketch of an unbound orbit of energy E > mc 2 and angular momentum L slightly 
larger than a/c in a 1/r potential. Because L is only slightly larger than a/c, the orbit makes 
loops around the scattering center. Such loops do not occur for an unbounded nonrelativistic 
orbit which is always a parabola or hyperbola. 

Fig. 3. Sketch of an unbound orbit of energy E > mc 2 and angular momentum L < a/c. 
One end of the orbit is at infinity and the other corresponds to a spiral into or out of the 
center of the potential. The nonrelativistic limit corresponds to a straight line orbit with 
L = 0. 

Fig. 4. Sketch of a bound orbit of energy < E < mc 2 and angular momentum L > a/c. 
The orbit has two turning points and, in the nonrelativistic limit, reduces to the familiar 
elliptical orbit of nonrelativistic mechanics. 

Fig. 5. Sketch of a bound orbit of energy E < mc 2 and angular momentum L < a/c. 
The orbit ends in spirals to and from the potential center. The nonrelativistic limit is a 
straight line orbit with L = 0. 
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